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Abstract 

The concept of system signature was introduced by Samaniego for systems 
whose components have i.i.d. hfetimes. We consider its extension to the 
continuous dependent case and give an exphcit expression for this extension 
as a difference of weighted means of the structure function values. We then 
derive a formula for the computation of the coefficients of these weighted 
means in the special case of independent continuous lifetimes. Finally, we 
interpret this extended concept of signature through a natural least squares 
approximation problem. 
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1. Introduction 

Consider an n-component semicoherent system. The design of such a 
system can be described through its structure function 0:{O,1}" {0,1}, 
which expresses at any time the state of the system in terms of the states of 
its components. Here "semicoherent" means that the structure function is 
nondecreasing in each variable and satisfies the boundary conditions 0(0) = 
and 0(1) = 1. 

Let Xi, . . . , Xn denote the component lifetimes and let Xi.n, . . . , Xn-.n be 
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the order statistics obtained by rearranging the variables Xi,. . . , in as- 
cending order of magnitude; that is, Xi.n ^ ••• ^ Xn-.n- 

The concept of signature was introduced in 1985 by Samaniego [8] for 
systems whose components have i.i.d. hfetimes as the n-tuple s = (si, . . . , s„) 
defined by 

Sk=Vl{T = Xu:n) (1) 

where T denotes the system hfetime. That is, Sk is the probabihty that the 
fcth component failure causes the system to fail. For general background, see 
Samaniego [9|. 

Under the i.i.d. assumption, the probability Sk can be computed as the 
ratio of 7?,^, the number of orderings for which the fcth component failure 
causes the system failure, to n!, the total number of possible orderings of the 
failure times. An alternative way to calculate s^, which does not require the 
exhaustive inspection of all the orderings, was found by Boland jl[ through 
the formula 

^^ = r^ E <^(x)-tIt e m (2) 

\n-k+l) x£{0,ir \n-k) X£{0,1}" 

|x|=n-fc+l |x|=n-A: 

where |x| = YJi=i ^i- 

Even though the i.i.d. assumption is quite reasonable in many applica- 
tions, especially when we want to compare different system designs, it is 
legitimate to investigate the probability ([1]) in the general case of depen- 
dent lifetimes and to search for formulas which extend ([2]) to this general 
framework. We observe that only a few results have been obtained in this 
direction, assuming for instance that the component lifetimes are exchange- 
able or independent and exponentially distributed (see, e.g., Samaniego |9|, 
§8.3] and Navarro et al. |6!]). 

Equation ([2]) shows that, in the i.i.d. case, the n-tuple s depends only 
on the system design. This is no longer true in the general dependent case, 
where the probability Pr(T = Xk-.n) may also depend on the joint c.d.f. of the 
component lifetimes, that is the function 

F(ti,...,t„) =Pr(Xi 

In this general setting, we shall denote the probability Pr(T = Xk-.n) by pk to 
emphasize that the ri-tuple p = (pi, . . . is not always a signature in the 
strict sense (i.e., a feature of the system design only). 
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In this paper, assuming only that the joint c.d.f. F is absolutely continu- 
ous, we provide a closed-form expression for = Pr(T = Xk.n) as a difference 
of two weighted arithmetic means of the structure function values whose 
weights depend only on F (Theorem [2]). We show that this expression re- 
duces to ([2]) as soon as the component lifetimes are exchangeable (Remark [2]). 
We also provide a useful expression for the weights (as a one- dimensional inte- 
gral) in the special case of independent lifetimes (Proposition [3]) and examine 
the particular case of independent Weibull lifetimes, which includes the ex- 
ponential model (Corollary H]). Finally, we show that the n-tuple p can be 
obtained from a symmetric approximation of the structure function in the 
sense of weighted least squares and we point out a formal analogy between 
this approximation and the concept of projected system recently introduced 
in Navarro et al. [tI]. 

Through the usual identification of the elements of {0, 1}" with the sub- 
sets of [n] = {l,...,^}, a pseudo-Boolean function /:{0,1}" ^ M can be 
equivalently described by a set function t>j:2["] R. We simply write 
Vf{S) = where I5 denotes the n-tuple whose ith coordinate is 1, 

if z € S*, and 0, otherwise. To avoid cumbersome notation, we henceforth 
use the same symbol to denote both a given pseudo-Boolean function and 
its underlying set function, thus writing /:{0,1}"' -> M or /:2["] M inter- 
changeably. 

The kth order statistic function oskm' {0, l}*^ {0, 1} is defined by oSfc:„(x) = 
Xk:n- We then have oSfc:n(x) - 1, ii \x.\ ^ n - k + 1, and 0, otherwise. As a 
matter of convenience, we also formally define oso:n = and os„+i:„ = 1. 



2. Explicit expressions 

Let F be the (absolutely continuous) joint c.d.f. of the component life- 
times Xi,. . . ,Xn. We define the associated relative quality function q:2^'^^ ->■ 
[0, 1] as 

q(S) = Pr ( max X,: < minX, ) 

with the convention that q{0) = qi[n]) = 1. That is, q{S) is the probability 
that the lifetime of every component in S is greater than the lifetime of every 
component in [n] \ S. Thus defined, q{S) is a measure of the overall quality 
of the components in S when compared with the components in [n] \ S. 
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Since the r.v.'s Xi,. . . ,Xn are continuous, we see that the function q can 
also be written as 

qiS)= E Pr(X.(i)<-<X.(„)) (3) 

creSn : {a(n-\S\ + l),...,a(n)}=S 

where &n denotes the group of permutations on [n]. 

We then observe that, for every k e [n], the values q{S) for \S\ = k sum 
up to one. In fact, by ([3]), we have 

E q{S)= E Pr(X.(i)<-<X.(„)) = l. (4) 

|5|=fc 

Remark 1. (a) The validity of (jl]) is especially transparent if one focuses 
on a particular simple case, say, when k = 1. It then simply says that 
there exists (with probability 1) a unique component j whose lifetime 
is maximum, that is, 

y Pr max X, < Xj) = 1. 

j^l Vje[n]\{j} / 

(b) If the variables Xi, . . . ,X„ are exchangeable, then the relative quality 
function q is symmetric. In this case, by (jl]), we have q{S) = 

(c) Equation shows that comparing q{S) with q{S') is relevant when- 
ever \S\ = \S'\. In general, according to (b) above, the relative quality 
could be better measured by the function q{S) = (|^|) q{S). 

We now give an expression for pk = Pt{T = Xk-.n) as a difference of two 
arithmetic means of the structure function values weighted by the relative 
quality function. We first present a lemma. 

Lemma 1. For every k e [n], we have 

FTiT^Xk..n)= ^W0(X). 

|xj=n-fc+l 

Proof. For every k e [n] and every a e we have 

Pr(r ^ X.(fc) I < - < = (l){{a{k),. . . , a{n)}). (5) 
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Indeed, assume ^ct(i) < ••• < X„(^n)- Using the path representation of the hfe 
function [Ij of the system, we have T ^ Xo-(fc) if and only if 

max minXj ^ Xfj(k\ (6) 

where Pi, ... , Pm denote the minimal path sets of the system. In turn, event 
([6]) is equivalent to 

"There is a minimal path set P such that Xi ^ ^^-(jt) for all i e P" 

Equivalently, "there is a minimal path set P £ {(j{k), . . . , a{n)}" . By mono- 
tonicity of 0, this event reduces to "0({(7(A;), . . . ,(T(n)}) = 1", which finally 
proves ([S]). 

By combining the law of total probability with we obtain 
Pr(T^X,,„)= E 0({^(A;),...,a(n)})Pr(X.(i)<-<X^(„)). 

Grouping the terms for which {a(k), . . . ,(j(n)} is a fixed set S and then 
summing over S, we obtain 

Pr(T^Xfc.„)= E E (/*(5')Pr(X,(i)<-<X,(„)). 

|5|=n-fc+l o-e6„:{o-(A:),...,(T(n)}=S' 

The result then follows from ([3]). □ 
Theorem 2. For every k e [n], w;e have 

p,= E 9(x)0(x)- E ?(x)</>(x). (7) 

|x[=n-fc+l |x|=n-fc 

Proo/. We have Pr(T = Xt.n) = Pr(T ^ Xfe:„) - Pr(T ^ We then 

conclude by Lemma [H □ 

Remark 2. (a) It is noteworthy that pk can be rewritten in the form 

Pk= E '"fc(x)0(x), 

xe{0,l}" 

with rfc(x) = q'(x) + 2 2;,fc:„ - Xk-i-.n)- This fact follows immedi- 

ately from the identity 

E g(x)0(x)xfc:n = E 9(x)</'(x). 

xe{0,l}" |x|jn-fc+l 
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(b) As expected, we observe from ([7]) that pk combines linearly two pieces 
of information: 

(i) the system design, which is encoded in the structure function 0, 
and 

(ii) the component lifetimes, which are encoded in the relative quality 
function q. 

(c) When the variables Xi, . . . ,X„ are exchangeable, by Remark W^b) we 
see that ([7]) reduces to ([2]). Thus in this case, p is the signature of the 



se, p 

system. This fact was previously observed in p. Lemma 1]. 

The following proposition gives a formula for q{S) as a one- dimensional 
integral in the special case of independent lifetimes. 

Proposition 3. For independent continuous lifetimes Xi, . . . ,Xn, each Xi 
having p.d.f. fi and c.d.f. Fi, with Fi(0) = 0, we have 

m n m n ii-m))dt. 

Proof. Denote the p.d.f. and c.d.f. of a r.v. X by fx and Fx, respectively. 
Recall that for two independent continuous r.v.'s X and Y, we have 

Pr(X<y)= [f fx{x)fy{y)dxdy= Fx{y)fY{y)dy. 

JJx<y J —oo 

By applying this formula to q{S), we obtain 

qiS) = f Pi( max Xi^t)— Ft (minXi^t) dt 

= r n m^ii-m-Fm]dt 

which immediately leads to the result. □ 

Corollary 4. For independent Weihull lifetimes, with Fi{t) = 1 - e^^'*'**)", we 
have 

for every S 0, where Xa{S) = T,ies^i- 
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Proof. By Proposition El we have 

ie[n]\S 

where the product can be expanded (by the generahzed binomial theorem) 



as 



We then have 

which immediately leads to the result. □ 

Remark 3. (a) Given a set function g:2["'] [0, 1], there is an additive set 
function Aq:2["'] ^ ]0, cx3[ satisfying (|8]) if and only if 

q{[n]\{i})>0 Vie[n] (9) 

and 

Indeed, if such an additive function exists, then by ([8]) we obtain 

«(N^W) = 1-^44^ = ^ (11) 

A„([nJ) A„([nJ) 

which leads to (|9]) and (fTOj) . Conversely, if (|9]) holds, then we may 
choose Xa{{i}) = q([n] \ {i}) and we then see that (fTOl) leads to (|8]). 
Thus ([9]) and (|TOl) provide necessary and sufficient conditions on a set 
function q:2^'^^ [0, 1] to be a relative quality function obtained from 
Weibull lifetimes (with a common shape parameter a). 
(b) Under the assumptions of CorollaryHl by f lTT]) the ratio AQ,({i})/AQ,([n]) 
is exactly the probability that Xi is the shortest lifetime. More gener- 
ally, the ratio Xa(S) / Xa{[n]) = T.ieslil''^] ^ {0) is the probability that 
the component having the shortest lifetime is in 5*. 
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(c) For S £ [n], the S'- difference of a function /: {0, 1}" ^ M is defined in- 
ductively by A^/ = / and A^/ = Ai^}A-5-{*}/ for i e S, with A{^}/(x) = 
/(x \ Xi = 1) - /(x I Xi = 0). It is then easy to see [3|, §2] that ([8]) can 
be rewritten as 

qiS) = (-ir'^l X^{S) (aN^^^^)(^). (12) 

Moreover, the {[n] \ 5) -difference in (fT2l) can be interpreted as the 
marginal interaction [3, §2] (associated with the function l/A^) among 
the components in [n] xS" conditioned to the presence of the components 
in S. 



3. Links with approximations of structure functions 

In Q the authors solved the problem of approximating a given pseudo- 
Boolean function /: {0, 1}" M by a symmetric one in the sense of weighted 
least squares. 

Specifically, given a weight function w:{0,l}" ]0,oo[, the best sym- 
metric approximation of a function f:{0,l}^ -> M is defined as the unique 
symmetric function /*:{0,1}"->M that minimizes the weighted squared dis- 
tance 

\\f-gr= E ^x)(/(x)-^(x))' (13) 

xe{0,l}" 

among all symmetric functions g: {0, 1}" -> M. 

The best symmetric approximation /* is actually the orthogonal projec- 
tion of /, with respect to the inner product 

{f,g)= ^(x)/(x)^?(x) 

xe{0,l}" 

onto the linear subspace of symmetric functions g-{0, 1}" M. In terms of 
the order statistic functions, this projection is given by 

n 

r = f{0) + Y.CkOSk:n (14) 

fc=l 

where 

c,= E «J(x)/(x)- E ^(x)/(x) (15) 

|x|=ri-fc+l |x|=7i-fc 
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and 

^i;(x) 
w(x) = — 

which shows that Ck is actually a difference of two expected values (see 
Since Ck is the coefficient of oSfc:„ in /*, it can be interpreted as a measure of 
the influence of the kth smallest variable on /. 

Now, consider an n-component semicoherent system defined by a struc- 
ture function (j) and an absolutely continuous joint c.d.f. F of the component 
lifetimes. Assume that the associated relative quality function q is strictly 
positive. Consider also the weighted distance (fT3ll with w - q and apply the 
approximation problem above to the structure function f - (j). By ([1]) we see 
that w = q = q. Theorem [2] then shows that the coefficient c^, as defined in 
f|T5|) . is precisely the probability pk = Pr(T = X^-.n)- 

Moreover, from f|T^ it follows that the best symmetric approximation (p* 
of (p (with respect to the weighted distance ffT^ with w = q) is given by 

n 

4>* = Y.PkOSk:n- (16) 
fc=l 

Under the i.i.d. assumption, (fT6!l reduces to 0* = ELi^fcOs^in, where Sk is 
given by ([2]). 

Interestingly enough, we also observe a formal analogy between the or- 
thogonal projection f|T6|) of and the concept of projected system, recently 
introduced in Navarro et al. Indeed, (j)* is a combination of the order 
statistic functions weighted by the n-tuple p while the projected system is 
the system which mixes the A;-out-of-r;, systems (0^ = osk-.n) with mixing dis- 
tribution p. 
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